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Kinetics of dislocations is studied by means of computer simulation during intensive plastic de- 
formation. The dynamical effect in the form of soliton-like wave of sharply disrupted interparticle 
bonds is observed. Along with it, micropores similar to steam bubbles at water boiling are formed. 
After some deformation the solid takes an ideal structure again owing to "emission of the bubbles". 
The dislocations in this state (with micropore) can form a grain boundary in the case of nano- 
structural materials. The nanoscopic object of rotational nature is observed at uniaxial intensive 
deformation. Rotation of some volume of sample which realized with phase transition scenario was 
observed by intensive shear deformation. 

PACS numbers: 62.20.Fe; 62.50.+p; 63 

Keywords: molecular dynamics, dislocation, intensive plastic deformation, soliton-like wave, phase transition 



I. INTRODUCTION 

A new methods of intensive plastic deformation for 
achievement of nano-structure in metals had been devel- 
oped §. It is common knowledge that the main vehicle 
of plastic deformation is dislocation. In recent papers 
new vehicle of plastic deformation had been investigated 
widely for metal glasses. It is so called shear transforma- 
tion zone (STZ) f|, §, @, @. On other hand, there is 
a large body of research of different dynamical objects in 
non-linear chains. They are moving breathers in differ- 
ent non- linear systems || and kinks in acoustic chains , 
M. The main property of these objects is spontaneous 
concentration energy in the large excitations at cost of 
suppress of weak ones. It is possible to assume that such 
concentration in three-dimensional case can lead to rear- 
rangement of internal structure of solid during intensive 
deformation. 



with Cartesian coordinates X^j, and X\ k , Y\ k . Indexes 
i, I numerate atoms in the lattice along the ^-direction, 
j, k do the same along the A"-direction. Let use a unit 
in which the constants are equal to A = 0.0025J c mJ? , 
B = lmj?, where J c and m c are conventional units of en- 
ergy and length respectively. We put the mass of a par- 
ticle to be M — 0.01fc<7 c , where kg c is conventional unit 
of mass. The binding energy of pair interaction is found 
to be E b = A/12B = 0.00020833J C and equilibrium dis- 
tance is ?'o = v^B = lm c . In case of real material it is 
necessary to perform the re-computation with account of 
real values of Eb, ro and the mass M. For example, in the 
case of copper E b = 0.5493 * 1CT 19 J, r = 2.66 * 10~ 10 ra, 
M = 1.0541 * 10~ 25 fc3. Then, lm c = 2.66 * 10" 10 m, 
lkg c = 1.054 * 10 - 23 kg, ls c = 0.532 * 10~ 13 s. The time 
unit have the same order as the period of small vibrations 
T = 2nm/U'lj lm . Time step for simulation is chosen to 
be At = 0.18s c . 



II. CONSTITUTIVE RELATIONS 



III. UNIAXIAL LOADING 



For the check of mentioned above assumption numeri- 
cal experiment on two-dimensional model was performed 
by means of MD-simulation. Model system of 18*18 par- 
ticles with Lennard-Jones potential is used: 
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where r^ft = v(AV, - X lk ) 2 ~ (Yy - Yi k ) 2 - is the dis- 
tance between the particles numbered by i,j and l,k 
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The first series of real time computer simulations is 
fulfilled for uniaxial loading of a crystallite sample. The 
system is placed on rigid platform constructed from the 
atoms identical with those in the sample (the lowermost 
line of atoms on the figures). The rigid atomic line 
(plane in a three-dimensional case) constructed from the 
same atoms moves downward with a constant velocity 
2.778 * 10~ 4 m c /s c (the uppermost line of atoms on fig. 
|l|). Initial atom configuration is hexagonal lattice with 
interatomic distance ro = lm c . 

In the first stage the usual uniaxial elastic deformation 
goes (Fig.^Ja). Further, two microscopic defects occur al- 
most simultaneously in the left upper corner of the lattice 
and in the right down one (Fig.|l|b). Their Burgers vec- 
tors have identical direction (at a 120° angle with hori- 
zontal axis) and opposite signs. These defects are similar 
to usual dislocations |)| . However each of them has two 
and only two inserted atomic lines. They can't have one 
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FIG. 1: Evolution of crystallite structure at uniaxial loading: a - ideal crystal structure at time step 27007, b - birth of 
dislocation pair at 29253, c - full formed bidislocation pair at 35 35019, d - lightning-like breakthrough between bidislocations 
at 35923, e - turned bidislocation pair at 36042, f - pressing out of the bidislocations from the sample at 36117, h - graphics of 
total potential energy of the system and kinetic and full ones during the turn of the bidislocations. Average potential energy per 
atom is equal to — 0.00132 J c . Yellow circles present the atoms with potential energy 0. 00002 J c < U < 0. 00004 J c , green circles 
- with 0.00004J C < U < O.OOOUc, blue circles - with 0.0001J C < U < 0.00015J C , red circles - 0.00015J C < U < 0.00038J C ,. 
Potential energy of the particles on free lateral surfaces is equal 0.0006J C . 
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inserted atomic line since atomic configuration doesn't 
permit it owing to very large energy of the defects. Two 
inserted lines and glide line are oriented at 60° angle each 
to another. Both of the inserted lines always lie on one 
and only one side of the glide line. It is pertinent to 
call a defect with two inserted atomic lines as bidisloca- 
tion. Along with two-dimensional lattice, such defects 
may be observed in three-dimensional case on the planes 
of closely packed layers (planes {1,1,1}). In general case 
of hexagonal lattice six type of bidislocation having Burg- 
ers vectors with the angles 0, it, ±7r/3, ±2-7r/3 may take 
place. In complex plane these vectors are expressed as 



±a ,a exp(±-i),a exp(±— i), (2) 

where ao - is equilibrium distance in the crystallite. 

In these terms, arisen bidislocations have following 
Burgers vectors: 

TT 277 

aQexp(——ij and aoexp(—i). (3) 

Under following loading the bidislocations move along 
their glide lines for as long as they reach the rigid lines 
of atoms (Fig. |l|c). Then, the bidislocations stop for 
some time because the rigid lines play the role of obsta- 
cles for further moving. Owing to this retardation, the 
bidislocationons form the pair with zeroth total Burgers 
vector. In this case the rotation of the region between 
the bidislocations takes place. This region has approxi- 



mately circle form. The bidislocations lie exactly on di- 
ametrically opposite sides of the circle boundary. In this 
nanoscopic region the orientation of atomic lines (planes) 
along Burgers vector changes continuously. Owing to this 
property, observed pair of bidislocations may be simplest 
two-dimensional model of displanation JhJ, JTT) . The 
displanations are known as defects of mesoscopic level. 
It is clear from presented experiment that simplest mod- 
ification of the displanations has nanoscopic nature. 

Stopping of the bidislocations creates difficulties for 
repairing of ideal crystal symmetry. As a consequence, 
local concentration of elastic energy around defects takes 
places. It is clear that for finishing of symmetry repairing 
gliding lines must turn to the direction along rigid lines 
of atoms. What scenario does the system choice for it? 

At the time of 600 time steps flashover-like puncture 
of the region between the bidislocations occurs (Fig. |l]d) . 
After it both the bidislocations are turned (Fig. [ljs) and 
pressed out from the lattice in the direction parallel to 
rigid lines (Fig. |l]f). Thereafter, the sample takes ideal 
crystal structure. Then, with continuing of loading the 
processes are qualitatively repeated. 

It is known that different transformations of disloca- 
tions can be regarded as chemical |l^] or, better still, as 
elementary particle reactions. From this standpoint, in 
the case of two-dimensional hexagonal lattice any bidislo- 
cation may be considered as an elementary particle with 
vector charge (2) or as three particles (like quarks) with 
different own polarization and with internal coordinate 
(as spin). Let us denote the set of "elementary particles" 
as 



nf = ±a , irf = (ao exp(— i), a exp( 



Then the transition with rotation of glide lines (fig. 
l.d-f) can be written in the form of reaction: 

7T+ +7T 2 - -> 7T+ +7Tf . (5) 

Therewith, the law of conservation of charge is fulfilled: 

aoexp(— i) + ao exp(— — -i) = ao — ao = const = 0. (6) 
o o 

The law of conservation of charge momentum is ful- 
filled since the distances between glide lines pre and post- 
reaction are equal (7 lines between them on the fig. |l|d,f). 
The law of conservation of energy is fulfilled since the 
system is Hamilton one. As own energy of each of two 
bidislocations pre and post-reaction doesn't change (in 
consequence of their identic structure) accumulated elas- 
tic energy almost completely turns into kinetic energy 
(relaxation of stress) . 



-i)), irf = (a exp(-i),a exp( — — i)). (4) 
I 

Regarding all the pictures presented in the fig. [I], 
one can note that bidislocation nucleus isn't localized in 
the point, but it is continuously smeared along the glide 
line in the region of several atomic sizes. The discrep- 
ancy of crystalline atomic lines begins with the distor- 
tion of them at bidislocation edges and increases gradu- 
ally toward center. Near a glide line the chain of atoms 
viewed from the side of more rarefied region of crystallite 
is stretched. Let call this atomic chain as "stretched" 
chain. The chain viewed from the opposite side of the 
glide line is termed as " compressed" chain. The atoms in 
the stretched chain are localized in energetically unfavor- 
able positions to ones in the compressed chain. However, 
the atoms are fixed in these positions owing to effective 
interaction with those in their own chains and with near- 
est atoms of the rest lattice. 

Let set as 1 equilibrium interparticle distance in the 
region far from the bidislocations. Then, in the com- 
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pressed chain interparticle distance is equal to 1.1 and 
to 0.95 in stretched chain. Moreover, all the simulations 
exhibit the next trend. In the chains parallel to com- 
pressed and stretched ones and viewed from the com- 
pressed chain interparticle distance takes its asymptotic 
value 1 rapidly, already in the first chain. Alternately, 
in parallel chains viewed from the stretched one several 
chains are markedly stretched. Thus, the region of com- 
pression is more localized then the region of rarefying. 
Owing to this, potential energy of any particle in the 
compressed chain is more large than one in asymptotic 
zone. This peculiarity is the main sign for structure de- 
fects of any nature and permits to separate bidisloca- 
tions immediately during computer simulation. Moving 
colored bidislocations against the background of the rest 
atoms is very nice picture. Average potential energy per 
atom in the compressed chain consists of 60-70 percent 
of one on free surface of the sample. Potential energy of 
whole the system in initial state of lattice is chosen as 
zeroth level. 



IV. "KNIFE" LOADING 

More interesting result is obtained in the second se- 
ries of computer simulation for cutting the sample by 
three-atom knife (fig. ||) . Three atoms only are left from 
moving upper rigid line of atoms. Other parameters of 
the experiment remain as in the Sec. III. 

In the first stage, in general, the picture qualitatively 
looks like to one obtained in the previous numerical ex- 
periments. Elastic stage of deformation isn't so interest. 
Further, the pair of bidislocations is created in the upper 
left and right corners of the crystallite. Their glide lines 
are turned by 60° each to another. The bidislocations 
occur in the region under three-atom indentor and move 
downward. Then, the bidislocations turn and they are 
pressed out from the sample. 

Interesting phenomenon is observed at such pressing 
out (Fig. ||). The defect displayed on the figure has two 
additional lines and Burgers vector as the bidislocation 
does, but its form sufficiently differs. The stretched chain 
has gap in the center of bidislocation nucleus. That is, 
binding among centered atoms is broken and micropore 
arises. 

Owing to the micropore, the bidislocation becomes 
more local and more sharp. Atomic lines in the regions 
separated by the normal to glide line arc inclined at a big 
angle each to other. One can regard the line which sepa- 
rates such regions as an element of a boundary between 
grains. 

During numerical experiment new bidislocations with 
micropores arise and disappeare again and again. Whole 
the picture is similar to formation of steam bubble at 
water boiling. Owing to this "emission of the bubbles", 
the solid take on ideal structure again after some defor- 
mation. 

On the other hand, the movement of a such object 



through the lattice is similar to kink-like excitation in 
acoustic chain. Wave-like disruption and repairing of 
bonds between atoms give rapid pressing out of defects 
from a sample. 

Thus, at intensive influence on metal it may release not 
only by well known shear deformation, but by wave-like 
(or excitation-like) bond disruption, which moves rapidly 
enough. One can assume that the micropore, namely, 
is the sign of new state of bidiislocation, which arises 
just during intensive deformation. In three-dimensional 
case the micropores can be spherical or elliptic bubbles. 
After stopping of deformation they disappear through 
the "condensation" into usual dislocations. Probably, it 
is difficult to detect the bidislocations with micropores in 
experiment as separate objects. Usual X-ray diffraction 
analysis isn't suitable for investigation of such states and 
structure re-arrangements for two reasons. By the first, 
such analysis is fulfilled much later after finish of intensive 
plastic deformation. By the second, the solid in such 
process can be in non-equilibrium amorphous state with 
large widening of X-ray spectrum line. For detection of 
such states and its changes, it is necessary to execute X- 
ray analysis just during intensive plastic deformation or, 
maybe, immediately after it. 



V. SHEAR LOADING 

A great amount of such objects occurs in the case of 
intensive deformation by shear scheme. This scheme is 
realized in the following way: the sample consisted from 
28*28 particles is placed in rigid two-dimensional box 
constructed from atoms identical with those of sample. 
The box is deformed by horizontal movement of upper 
and down rigid lines with constant velocity. 

The created dislocations with micropores form bound- 
ary separated regions with different orientation of crys- 
talline atomic lines (Figjjj). Such boundary can play 
important role for forming of nano-structural state by 
means of the methods of large shear deformation. 

At the beginning stage of deformation single bidislo- 
cation is borne in the left down corner of the sample 
(Fig. ||a). It has the charge a exp(— 2ni/3) and moves 
along glide line toward the upper boundary. At the fol- 
lowing stage, the second bidislocation with the charge 
ao exp(— 7rz/3) comes into being and moves along glide 
line toward the upper boundary too (Fig. ||b). The glide 
lines of both bidislocations are crossed at the same point 
exactly on the upper boundary. With the lapse of time 
the first bidislocation takes the form of microcrack (Fig. 
||c) and preserves this form for a time (Fig. ||d). In the 
following lapse the first bidislocation begin to take the 
form of sharp bend (Fig. Us). Upon it, pair of bidislo- 
cations arises instead the first one. One of them has the 
charge — ao and another has the charge ao exp(— 7rz/3). 
Their total charge is equal to ao exp(— 2-7rz/3), i.e. it is 
equal to the charge of the initial first bidislocation. Such 
transformation can be written as decay reaction in the 




FIG. 2: Evolution of crystallite structure at cutting by three atom knife: a - ideal crystal structure in beginning stage, b,c - 
formation of usual bidislocation, d, f-h - formation of bidislocation with micropore, e - renewal ideal crystal structure after first 
"bubble". 
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FIG. 3: Evolution of crystallite structure with shear deformation: a - birth of first bidislocation, b-d - moving of two defects, e 
- bidislocation decay with birth of two new bidislocation, f - system after bidislocation decay. 
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form: 

+ energy — > ir^ + ir^ . (7) 

where the second term in the left part of equation is the 
energy of stress field which causes the birth of additional 
" particle" . The rest of stress field energy transforms into 
thermal energy. As the result, almost full relaxation of 
elastic stress take place. 

The dislocations with micropores form boundary sep- 
arated regions with sharply different orientation of crys- 
talline atomic lines (fig. |^). Such boundary can play 
important role for forming of nano-structural state by 
the methods of large shear deformation. Single, but very 
interesting result had been fixed in computer experiment 
when rotation of certain volume by scenario of phase 
transition takes place. At some deformation the bound- 
ary is formed as a chain of bidislocations with microp- 
ores. At following loading it moves steadily until escape 
from the sample. In this case the boundary separates 
phases of identical nature, but with different orientation 



of crystalline atomic lines only. 

With my point of view, the defects described above 
have similar nature to them considered in ref. || - ||. 
Our "STZs" differ from § - § ones because the former 
are observed in ideal crystallite, but no in amorphous 
metallic glass. 



VI. SUMMARY 

Thus, quark- like behavior of bidislocations is observed 
in different numerical experiments. It is established that 
the bidislocations may have at least two different struc- 
tural states. The first of them has whole stretched chains, 
the second has micropores. The bidislocations with mi- 
cropores, as a rule, have pentagon form. The states with 
crack-like form of bidislocations are possible too (Fig. 
||c,d). The transient bent form of microcrack arises at 
the decay of dislocation (Fig. |^e) . The superhigh-speed 
soliton-like movement of bidislocations with micropore is 
noted. 
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FIG. 4: Phase-like transition with moving boundary consisting from bubble-bidislocations having pentagon form. 



